Euler Integration of Gaussian Random 
Fields and Persistent Homology 

Omer Bobrowski* >t >+ , Matthew Strom Borman* 

Omer Bobrowski 
Department of Electrical Engineering 
Technion, Israel Institute of Technology 
Haifa, 32000 
Israel 

e-mail: boberOtx.technion.ac.il 



Matthew Strom Borman 
Department of Mathematics 
University of Chicago 
Chicago, Illinois 60637 
USA 

e-mail: borman@math . uchicago . edu 

Abstract: In this paper we extend the notion of the Euler characteristic 
to persistent homology and give the relationship between the Euler inte- 
gral of a function and the Euler characteristic of the function's persistent 
homology. We then proceed to compute the expected Euler integral of a 
Gaussian random field using the Gaussian kinematic formula and obtain a 
simple closed form expression. This results in the first explicitly computable 
mean of a quantitative descriptor for the persistent homology of a Gaussian 
random field. 

AMS 2000 subject classifications: 60G15, 55N35. 

Keywords and phrases: Persistent homology, barcodes, Betti numbers, 
Euler characteristic, random fields, Gaussian processes, Gaussian kinematic 
formula. . 



1. Introduction 

The Euler characteristic, denoted by x(X), is an invariant of nice topological 
spaces X, which is additive in the sense that there is inclusion-exclusion: 

X (A UB) = x(A) + x(B) - X (A D B). 

Since this is the heart of a measure, one can make sense of integration with 
respect to the Euler characteristic. This new integration theory has many of 
the properties of Lebesgue integration, such as a Fubini theorem and integral 
transforms (e.g. Fourier and Radon transforms [14]). Aside from the theoretical 
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aspects, the Euler integral has recently started to give rise to applications. For 
example, [4] presents a target enumeration scheme which employs the Euler 
integral to perform data aggregation in sensor networks. 

One major drawback of the original Euler integral is that it was only defined 
for integer valued functions. In [3], two extensions were given, called the upper 
and lower Euler integrals, which are defined for a large class of real valued func- 
tions. In this paper we will use the upper Euler integral, but everything we do 
has a lower integral analogue. This integral operator, while having some draw- 
backs (particularly a lack of linearity), possesses a compelling Morse theoretic 
interpretation (see [3]). In this paper we extend this Morse theoretic interpreta- 
tion to general tame functions. Additional information on the stratified Morse 
theory used in this paper can be found in [9], For algebraic topology prelimi- 
naries on the Euler characteristic and homology see [10]. 

We will use recent developments in the theory of Gaussian random fields 
to compute the expected value of the upper Euler integral of Gaussian and 
Gaussian related random fields. Random fields are stochastic processes defined 
over a topological space X. Gaussian fields are characterized by having finite 
dimensional distributions which are multivariate normal. A field generated by 
applying a function to a vector valued Gaussian field is called a Gaussian related 
field. Let / be a fc-dimensional, Gaussian field on a nice space M of dimension 
d, G : R fe —> R a piecewise C 2 function, g = Go/. Assume that / has zero 
mean and constant unit variance. Then, under mild regularity conditions, our 
most general result states that 




where D u = G _1 (— oo, u], E {g} := E {g(t)} (for any t e M), and the Cj,Mj are 
geometric characteristics known as Lipschitz-Killing curvatures and Gaussian 
Minkowski functionals, respectively. 

There is a strong link between the Euler integral and persistent homology. 
Briefly, the persistent homology of a real valued function / tracks changes in 
the homology of sublevel sets oo,it]. As the level u changes, new homol- 

ogy elements (i.e. "holes") are born and others die. Persistent homology keeps a 
record of this birth-death process. The theory of persistent homology was cre- 
ated as a way to describe the topology of data, but to date most results about 
persistent homology have been algebraic in nature as opposed to probabilistic 
or statistical. For the theory of persistent homology to develop into a powerful 
applied tool there will need to be theorems about persistent homology in the 
random setting. A first step in this direction was [5], where the expected bar- 
code of the persistent homology for points sampled from a circle was computed. 
We will relate persistent homology to the Euler integral via a parameter we 
call the Euler Characteristic of the Persistent Homology. Thus, we obtain the 
first known general tool that allows one to make probabilistic statements about 
the persistent homology generated by sublevel sets. Furthermore, our compu- 
tation of the expected Euler integral for Gaussian random fields provides the 
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first quantitative measure of the persistent homology of random fields in the 
multi-dimensional setting. 

Another significant and surprising corollary of the formula (1.1) is related to 
the signed sum of critical values of Gaussian random fields. If / : M — > R is a 
Gaussian random field and M is a closed manifold we prove that 

J (-ir (p) f( P ) 

where C(/) is the set of critical points of the field /, and fi{jp) is the Morse 
index of / at p. In other words, the expected signed sum of critical values of a 
Gaussian field does not scale according to the volume of the space as one might 
expect, but rather according to a one dimensional measure of the space. 

Finally, using (1.1) we can extend the application discussed in [4] to the 
situation where the sensor field is contaminated by additive Gaussian noise. In 
this particular case the integral remains additive, which enables us to suggest a 
noise reduction scheme. 



A(M) 
V2^ ' 



2. The Euler Integral 



The Euler characteristic is an additive operator on compact sets. Therefore, it 
is tempting to consider x as a measure and integrate with respect to it. The 
main problem in doing so is that x is only finitely additive. At first (see [14]), 
integration with respect to the Euler characteristic was defined for a small set 
of functions called constructible functions defined by 



CF{X)=\ f(x) = Y,a k l Ak (x) 



k=l 



<2fc G Id, A k disjoint tame subsets of X 



where 'tame' means having a finite Euler characteristic. For this set of functions 
we can define the Euler integral analogously to the Lebesgue integral. Let fix) = 
Sfe=i a k~&-A k {x), and define 



« n 

Jx 



k=l 

This integral has many nice properties, similarly to those of the Lebesgue inte- 
gral, such as linearity and a version of the Fubini theorem (see [14]). However, 
as mentioned above, the Euler characteristic is not countably additive, and 
therefore we cannot continue from here by approximating other functions using 
functions in CF(X). 

In [3] two extensions were suggested for the Euler integral of real valued 
functions using the notion of a definable function over an O-minimal structure 
(see [3] and references therein for more background). Let [^J and \x] be the floor 
and ceiling values of x respectively. In the O-minimal language, if / : X — > R is 
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a definable function on a definable space X, then an important property is that 
both [/J and [/] are constructive functions and hence have well defined Euler 
integrals. This leads to the following Riemann-sum like definition: 

Definition 2.1 ([3]). The lower Euler integral of a definable function f ; X — > 1R 
on a definable space X is defined by 

f[d X \ = lim n- 1 f [nf\d X , (2.1) 
x n ^°° Jx 

and the upper Euler integral is defined by 

f f\dx\ = lim n- 1 f \nftd X . (2.2) 
Jx n ^°° Jx 

These two extensions coincide with the original Euler integral for constructible 
functions. For other functions they might be completely different. 

Unfortunately, it is not clear if Gaussian random fields / : X — > K can be 
made to fit inside an O-minimal setting. Therefore in this work we will use the 
following simplified definition of a tame function. 

Definition 2.2. A continuous function f : X — » M on a compact topological 
space X with a finite Euler characteristic is tame if the homotopy types of 
/ (— oo, u] and f [u, oo) change only finitely many times as u varies overM. 
and the Euler characteristic of each set is always finite. 

We can now take our definition of the lower and upper Euler integrals to be 
as follows. 

Definition 2.3. If f : X — > M is a tame function, then the lower and upper 
Euler integrals are defined by 



f[dx\ = / lx(f > u) - X (f < -«)] du (2.3) 

rOO 

f\dx\ = / [X(f > u) - X (f < -«)] du, (2.4) 

Ju=0 

where x(f > u) := x (f^ 1 [u, oo)) , x(f < u ) ■= x( x ) - x(f > u ), etc. 

These equations appear as Proposition 2 in [3] , where they were proven using 
Definition 2.1. Thus, Definition 2.1 and Definition 2.3 agree when applied to a 
function that is both definable in the O-minimal sense and tame in the sense 
of Definition 2.2. Rather than viewing (2.3) and (2.4) as consequences of the 
O-minimal definitions (2.1) and (2.2), we will take (2.3) and (2.4) to be the 
definition of the Euler integrals for functions that are tame in the sense of 
Definition 2.2. 



2.1. The Euler Integral and Morse Theory 



In [3] the Euler integral was given a stratified Morse theory interpretation. A 
corollary of this approach was that if / : M — > R is a Morse function, and M is 
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a closed manifold, then 

/ f\dx\ = E ("^/(p), (2-5) 
M pec(/) 

where C(/) is the set of critical points of / and n(p) is the index of p as a critical 
point. In our language of tame functions, we have the following proposition. 

Proposition 2.4. Let f : X — > R be a tame function and let CV(/) be the set 
of values where the homotopy type of f~ l (— oo, u] changes (the critical values of 
f). Then 



f f\dx\= E A x(/.")". 
x vecv(f) 

where A x (f, v) = xif < v + e) — xif < v — e), for sufficiently small e, is the 
change in the Euler characteristic o// _1 (— oo,u] as u passes through the critical 
value v. 

Proof. Label the critical values CV(/) = {v\, . . . ,v n } in increasing order such 
that vi < ■ ■ ■ < Vi < < Vi+\ < ■ ■ ■ < v n . If Vk < u < vu+i, then via a 
telescoping sum 

X (f<u) = A x (f,v 1 ) + --- + A x (f ) v k ), 

Xif > «) = X(M) - X {f <u) = A x (f, v k+1 ) + ■ ■ ■ + A x (f, v n ). 
Therefore for u E [0, oo) and u ^ =t^_?, 

n 

Xif >u)= E A xif' v j) 1 [0,i; 3 ]H and 
j=i+i 

i 

Xif < -«) = E A xif,vj) Ip.-^lW- 
i=i 

Thus 

/ f\dx\= ixif >u)- xif <~u))du = Y. v ^xif^), 

JM Ju=0 j =1 

as desired. □ 

This recovers the Morse theoretic viewpoint, since if / : M —> R is a Morse 
function then Morse theory says that the Euler characteristic changes by the 
addition of (— 1) as /(— oo, u] passes through a critical point of index k. 

Corollary 2.5. Let f : X — > R be tame, satisfying the conditions of Proposi- 
tion 2.4- If a is not a critical value of f , then 

E A x (f,v)v= f G a {f)\dx\+axif<a)-ax(X) 

»ecvc/) X 

where G a (x) = min(a;,a). 
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Proof. If v\ < ■ ■ ■ < V n are the critical values of / and a is such that < a < 
Wfc+i, then G a (f) has critical values at v± < ■ ■ ■ < < a. By Proposition 2.4, 

/ G a (f)\dx\ =T / v j A x (G a (f),v j ) + aA x (G a (f),a) 

k 

= Y, vj A x (/, vj) + a( X (M) - X (f < a)). 

3=1 

This gives us the desired result. □ 

3. Gaussian Random Fields and the Gaussian Kinematic Formula 

There has been extensive effort over the past few years to study the sample 
paths of smooth random fields from a general Riemannian manifold M to R fc . 
Specific examples in which this approach has had practical importance occur 
when M is a 3-dimensional brain or a 2-dimensional cortical surface. The basic 
(random) geometrical objects studied were the excursion sets of the random 
fields, namely f~ 1 (D), for nice subsets D of M. k , and the tools for quantifying 
these sets were those of differential geometry. The theory of this subject has 
developed rapidly over the past few years (see [2, 11, 12]). One of its most 
powerful results is an explicit expression for the mean value of all Lipschitz- 
Killing curvatures (among them is the Euler characteristic) of excursion sets for 
centered (i.e. E {/(<)} = 0), constant variance, C 2 , Gaussian random fields. The 
result presented in [2] links random field theory with integral and differential 
geometry, and leads to approximations of important objects in probability and 
statistics, such as the exceedence probabilities P(sup M /(i) > u) [1, 13]. The 
main theorem in [2] is called the Gaussian kinematic formula (GKF). We state 
it here without getting fully involved with the details. 

Theorem 3.1 (The GKF, [2, Theorem 15.9.4]). Let M and D C M fe be regular 
stratified spaces with M compact and D closed. Let f — (/i, . . . , /&) : M — > M. k 
be a k-dimensional Gaussian field, with iid components all having zero mean, 
unit variance and such that with probability one fj is a stratified Morse function. 
Then 

dim ht—i r. 

E{A(/- 1 ( J D))}= £ 1 

3=0 L 

where: 

• Ci{-) is the i-th Lipschitz- Killing curvature of X (with respect to the metric 
defined in (3.2) below), 

• Mi(-) is the i-th Gaussian Minkowski functional, 

• [fc] = (fc) oi m" ' an d UJn * s ^ e v °l ume °f th e n- dimensional unit ball. 



(2n)-^ 2 C l+J (M)M j (D), 
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The definitions of regular stratified spaces, the Lipschitz-Killing curvatures, 
and the Gaussian Minkowski functionals are given in [2, Definition 9.2.2, Defi- 
nition 10.7.2, Corollary 10.9.6]. While the precise definition is involved, exam- 
ples of regular stratified spaces are closed manifolds, compact manifolds with 
boundary, and products of regular stratified spaces. Conditions are given in [2, 
Corollary 11.3.5] on the covariance function of /, such that sample paths are 
stratified Morse functions with probability one (and so in particular are tame 
functions). See [2] as well for a reference for most of the following facts. 

Lipschitz-Killing curvatures are geometric objects that depend on the choice 
of a Riemannian metric on M, such that Ck(M) is a measure of the fc-dimensional 
'size' of M. This means that if we scale the metric by A, then Ck(M) scales by 
A fe . For a large class of spaces, which include smooth manifolds and convex com- 
pact regions, if X C R™ is given the Euclidean metric, then the following tube 
formula holds for sufficiently small r 

n 

n(B(X,r)) =Y,"j£n-j(Xy , (3-1) 

where fj,(B(X,r)) is the Lebesgue measure of the tube of radius r about X. It 
turns out that the Cj(X) are independent of how X is isometrically embedded 
into R™ (and of n). For a d-dimensional space, Cd(X) is its Riemannian volume 
and Cq(X) is always its Euler characteristic. Here are few examples, if X C M 2 
is convex and compact then 

C (X) = 1, d(X) = (perimeter of X)/2, £ 2 (X) = area(X). 

Now suppose that M is a closed d-dimensional manifold. If M is odd dimensional 
then the even Ck vanish, while if M is even dimensional then the odd Ck vanish. 
When the parities match, Cd-j(M) for even j is given by an integral of an 
expression involving the curvature tensor. In the case that M has constant 
sectional curvature re, then 

£ d _ J .(M) = (4,)-^ (d _.^. /2) , ^Vol(M). 

The Gaussian Minkowski functionals (A4i(D)) are also defined using a tube 
formula similar to the one in (3.1), but using the Gaussian measure on M fe 
instead of Lebesgue measure. 

If / : M — >■ R is a mean zero Gaussian random field, then it is completely 
determined by its covariance function C : M x M — > R, where C(s, t) = 
E{/(s)/(<)}, and we can use this to define a metric g on M by 

g t (X t ,Y t ) =E{X t (f)Y t (f)} = X s Y t C(s,t) \ s=u (3.2) 

where X and Y are vector fields on M defined near t G M. It is with respect 
to this metric that the Lipschitz-Killing curvatures are computed in the GKF, 
and also the metric with respect to which we require M to be bounded. It turns 
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out that any Riemannian metric can be realized as coming from a variance one 
Gaussian random field. 

We will be particularly interested in the special case of the GKF for which 
i = 0. Since £q is just the Euler characteristic Xi the GKF implies that 

dim M 

E{x(/" 1 ( J D))}= £ (2n)-^ 2 C j (M)M j (D). 



4. The Euler Integral of Gaussian Random Field 

Let M be a stratified space and let g : M — > M be a Gaussian or Gaussian 
related random field. We are interested in computing the expected value of the 
Euler integral of the field g over M . While we focus on the upper Euler integral, 
everything we do has a lower Euler integral analogue. The following result is a 
corollary of the GKF and Proposition 2.3. 

Theorem 4.1. Let M be a compact d- dimensional stratified space, and let 
f : M — >• M. k be a k-dimensional Gaussian random field satisfying the GKF 
conditions. For a piecewise C 2 function G : R k -> K, let g = G o f . Setting 
D u = G _1 (— oo,u], we have 

e{ / g\dx\)^x{M)^{g}-Y,{^r ]/2 ^(M) [ Mj(D u )du (4.1) 

IJM ) ■ 1 JR 

where E{.g} = E{p(t)} (^(t) /ias a constant mean). 

The difficulty in evaluating the expression above lies in computing the Minkowski 
functionals Mj(D u ). In Sections 4.1 and 4.2 we present a few cases where they 
have been computed, which allows us to simplify (4.1). 

Proof. By Proposition 2.3 

/ sT^xl = / (xig >u)- xig < -u)) du 

Jm Jo 

= / ixiM) - xig < u )) du - / xig < u )du. 

JO J-oc 

Therefore, 

f ' g\dx\)= [°° (xiM)-E{ x (g<u)})du 





E{ x (g<u)}du. 

-oo 

Replacing D with D u in the GKF (Theorem 3.1) yields 

d 

E{x(.9 < u)} = E{x(r 1 (D u ))} =J2W j/ %iM)MjiD u ). 



(4.2) 
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Substituting this formula into (4.2) yields, 



/ g\d x -\\ = Y,W j/2 P^( M ^ 
Jm J J=0 



where 

- $ R M 3 {D u )du j>0, 
/ °° (1 - M (D U )) du - /? b (M (D u )du) j = 0. 



P.! 



The expression for po can be further simplified. Let X be a standard normal 
variable and Y — G(X), then 



Therefore, 



Mo (A,) = 7fc(AJ = P e A) = P (Y < u) . 



p = I {l-¥(Y <u))du- I P{Y<u)du 



(i 



' (Y > u)du— j F(Y <u)du 



o 



= E{r}. 

Since for all t, f(t) — Af (0,1), we can replace Y with G(f(t)) = g(t). Finally, 
recalling that £ = \ completes the proof. □ 



4-1. Real Valued Fields 

For real valued fields we can improve Theorem 4.1 by computing the terms 
■Mj(D u ) that appear in (4.1). First, we need to recall some facts about the 
family of Hermite polynomials. For n > 0, the n-th Hcrmite polynomial is 
defined as 

in 

H n {x) = {-lT V (x)^ — V {x), 

where (f(x) = (2n)^ 1 / 2 e~ x2 ^ 2 is the density of the standard Gaussian distri- 
bution. This family of polynomials is orthogonal under the inner product on 
functions /, g : K M 

(f,9) = / f{x)g{x)ip{x)dx. 
Jr 

A useful convention is 



/>oo 

H—i(x) = f(x)^ 1 / (p(u)du. 

J X 
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Theorem 4.2. Let M be a compact d- dimensional stratified space, and let f : 
M — > R be a real valued Gaussian random field satisfying the GKF conditions. 
Let G : R — > R be piecewise C 2 and g = G o f. Then 



f 9\dx\ } = X(M)E{ 9 } + J2 (-1) J '^(M) 
Jm J j=1 



(g 3 --i,( S ign(GQ)J^) 
(2tt)J/2 



In the case that the function G is strictly monotone, this can be simplified 
by using the fact that sign(G') is constant and then integrating by parts. 

Corollary 4.3. Let f be as in Theorem 4-. 2, and G be a strictly increasing 
function. Then 

'Hj,G) 



E {/ M ^l}=g™M)ig )j/2 



If G is strictly decreasing then, 

d 



-{/ M *l} = £^(M)ig 



Hj,G) 



j/2 

j=0 

To prove Theorem 4.2 we will need the following calculus lemma, which is a 
special case of Federer's coarea formula. 

Lemma 4.4. Let h : R — > R be an integrable function and let G : R — » R 

&e a piecewise differentiable continuous function that is nondifjerentiable on a 
discrete set. Then 



h(x) \G'(x)\dx 



[ I £ h(x)\dt. 



Proof. (Theorem 4.2) By Theorem 4.1, it suffices to show that 

M 3 (D u )du = (-iy (Hj-u (sign(G')) J G') , (4.3) 



for j > 1, where D u = G 1 (— oo,u]. Since G is continuous, we can write the 
inverse image of (— oo, u] as a disjoint union of closed intervals 

Du = [J[ai,bi] 

i 

where we allow one aj to be — oo and one bi to be oo. Note that for all the finite 
values we have G(ai) = Gipi) = u, G'(a,) < and G'ipi) > 0. 
For small enough p we have 

Tube (D u , p) = \j[a,i - p, h + p\. 
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Therefore 

7fc (Tube (D u , p)) = + P) - - p)) , (4.4) 

i 

where <J>(x) = (p(u)du. The Taylor expansion of <&(x + p) in p is 

OO j 

<f(x + p) = ${x) + ^(-ly-'H^^xMx), (4.5) 

3=1 J ' 

so in particular JAj (—00, x] = (— l)^ 1 Hj_i(x)ip(x). Therefore we conclude that 
for j > 1 

Mj {D u ) = {{-^Hi-iibiMh) + H^iaiMai)) . (4.6) 

i 

Note that if bi = 00 (or a; = —00) its contribution to the volume of the tube in 
(4.4) is independent of p (1 or respectively). Thus, it will affect only Ado and 
we can assume all the and bi in (4.6) are finite and hence Ui{°ij bi} — G _1 (w). 
If j is odd, then from (4.6) we have that 

Mj (D u ) = Y {Hj-i(k)(p(bi) + Hj-iiciiWai)) = ^ H^x{x)y{x). 

If j is even, then 

Mj (D u ) = {-Hj-i(bi)<p{h) + Hj-tiaiMai)) 

i 

For the case that j is odd, apply Lemma 4.4 to get 
Mh 



(D u )du= / 2J Hj-x(x)(p(x) \ du 



= J Hj-i(x)ip(x)\G'(x)\dx 
= {Hj^,\G'\) 

= (-ir 1 (^_ 1 ,(sign(G')) J G')- 

So we have proved (4.3), when j is odd. If j is even, a similar calculation gives 
the desired result. □ 



4-2. Vector Valued Fields 

When / is a vector valued Gaussian field, it can be difficult to evaluate the 
Minkowski functionals Mj. In two cases though, it is possible to compute the 
mean Euler integral. 
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4-2.1. The x 2 ease 

Let M be a compact <i-dimensional manifold. A x 2 field with k degrees of free- 
dom is of the form g = G o /, where / = (fx, ■ ■ ■ , fk) '■ M —> K fc is a Gaussian 
random field with iid components, and G(x\, . . . , Xk) — Y)j—-\ xf. 

Theorem 4.5. The mean Euler integral for a x 2 field with k degrees of freedom, 
with k > d, is given by 



E 



/ g \dx] ) = k C (M) - A A (M) + - A (M) 

J M J V 71 " M9) K 



Proof. First note that in this case, Mj (D u ) = Mj (G _1 (— 00, u\) — when 
u < since G is nonnegative. In [2, Section 15.10.2] it is shown that for k > d 
and j > 1 



M (D u ) 



d J 1 Pk(x) 



dxi- 1 



x k-i e -x 2 /2 



Therefore, 



Mj(D u )du = 



d^pkix) 



where Vk(%) - , . , 

K ' r(fc/2)2( fe - 2 )/ 2 

d^Pkit) 



du = 2 

V^T ■'<> 



10 dxi- 1 

Computing for j = 1, j = 2, d > j > 3, we have that 



dP~ 



tdt. 



> poo p/ fe+1 \ 

Mi{D u )du = 2 / p k (t)tdt = 2V2 -^fr. 

Jo 1(2) 

/•OO />OC 

/ M 2 (D u )du = 2 p' k (t)tdt = 2, 
Jo Jo 

and integration by parts yields 



A1./U>.,)./» = 2(^^t 



Finally, noting that E {</} = fc completes the proof. 



□ 



^.2.2. The F case 

Let M be a compact e?-dimensional manifold and let / : M —> M. n+m be a vector 
valued Gaussian field with iid components, 



G(x) 



Em 
i=l ; 



En 2 ' 
j=l X m+i 

and g = G o f. In this case, it is proved in [2, Theorem 15.10.3] that for j > 1 
(G _1 [u, 00)) = + ™ 



E E * 

1=0 i=0 



m,n,j,l,i 



(mu 



\ n 
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for a set of constants C m ,nj,lj.- 

Using basic calculus we can show that for n > j : + 2 and for all m, the integral 
J °° A4j (G _1 [w, oo)) du converges. This can be used to compute the expected 
lower Euler integral J M g \dx\ rather than the expected upper integral that we 
have computed so far. Thus, we can conclude that for n > d + 2 the expected 
lower Euler integral is finite. For each n, m it is possible to compute the exact 
value, but no general formula is known. In order to compute the upper Euler 
integral, we need to compute Aij(G~ 1 (— oo,u}). We note that this is feasible, 
but technically too complicated to be pursued here. 



5. Weighted Sum of Critical Values 



Taking G(x) = Hi(x) — x in Theorem 4.2 and using Proposition 2.4 yields the 
following compact formula. 

Theorem 5.1. Let f : M — > R be a Gaussian random field satisfying the 
conditions of the GKF, then 

e\[ f\d x ]\=E{ £ A X (/, V )4=-^P (5-1) 
Um J Ucv(/) J V2n 

where CV(/) is set of critical values of f and A x (/, v) is the change in the Euler 
characteristic of oo, u] as u passes through v from below. In the case that 

M is a closed manifold, then 

4 e (-ir ( ^(4=-^P- ( 5 - 2 ) 

Uc (/) J ^ 

where C(/) is the set of critical points of f , and /J>(p) := fi(p, f) is the Morse 
index of the critical point p. 

In the case that M is a closed even dimensional manifold, C\(M) = so (5.2) 
states that 

J j2 (-ir {p) f(p)\=o. 

UeC(/) J 
This fact has the following alternative proof, namely: 



.pec(/) J Uec(-/) 

| E M) m(p,/) /(p) 
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The first equality holds because / and — / have the same law. The second 
equality holds deterministically, using that M is even dimensional, because the 
negative of a Morse function is still a Morse function and critical points of index 
\i become critical points of index d — [i. 

The thing to note about Theorem 5.1 is that the expected value of a weighted 
sum of the critical values scales like £\(M), a 1-dimensional measure of M and 
not the volume Cd(M), as one might have expected. Consider the following 
example: Let / : R d —> R be a Gaussian random field with covariance function 

C : M. d x M. d — » M given by C(x, y) = e~~ This covariance function induces 

the Euclidian metric on M. d and Theorem 5.1 implies that 

[o,L] d J V 2ir V Zir 

In comparison to Theorem 5.1, letting G(x) — x d and using Theorem 4.2 
we get that E { J M f d \d\ \ } depends on the volume Cd{M) (as well the other 
measures). So while in general the behavior of the critical points and the critical 
values depends on the volume, when one takes the weighted sum of the critical 
values a lot of cancellation occurs and the result only depends on a 1-dimensional 
measure. 

The result in Theorem 5.1 can be generalized to the case where we consider 
only critical values below some level a. Observe that taking a — > oo in the 
theorem below recovers the result in Theorem 5.1. 

Theorem 5.2. Let f : M M. be a Gaussian random field satisfying the 
conditions of the GKF. Then 



El a x(/>«W =-<p(a)£o(M) 

k c <r J (5.3) 
d 

3=1 

In the case that M is a closed manifold, then the left hand side above can be 
replaced with E |X) pe c(/): /( P )<a( _1 

Proof. We need to investigate the Euler integral of G a (f), where recall that 
G a (x) = xl(_ 00i0 )(a;) + aTL[ at00 )(x) is the cutoff function from Corollary 2.5. 
According to Corollary 2.5, 



eJ m/^)4=4/ 

.60V(« UA 
K v<a ) 



Ga(f)\dx\ }+aE{ X (f <a)}-a X (M). 

M ) 

The first term on the right hand side is computed in Lemma 5.3 and the second 
term is given by the GKF (Theorem 3.1). □ 



Bobrowski and Borman/Euler Integration of Gaussian Random Fields 



15 



Lemma 5.3. Let f : M — > M be a Gaussian random field, satisfying the GKF 
conditions. Then, 

e( / G a (f)\d X ]\ = X (M)(a-a*(a)-tp(a)) 



A I 



d 

-^(a)^(2 7 r)-^ 2 £ J -(M)/f J -_ 2 (-a) 
i=i 

Proof. We will apply the result from Theorem 4.1, so we need to compute 
E {G a (/)} and f R M j (G- 1 (-oo,u])du. Note that 



G- 1 (- W , w ] = 



(— oo, u] u < a, 
M u > a. 



For j > 1, by (4.5) we know that Mj((-oo, u]) = 1 H -i {u)ip{u) and 

jVfj(E) = 0. Therefore if j > 1, then using that H k (-x) = (-l) k H k (x) 



M 3 (D a u )du= / (-ly^Hj-^tpi^du 

I J~oc 

/oo 
Hj-i(u)<p(u) du 
-a 

where the last transition is due to integration by parts. We also have that 

/a p oo 

X(p(x)dx + / aip(x)dx = a — a$(a) — <p(a). 
-oo J a 

Thus, by Theorem 4.1, we are done. □ 



6. The Persistent Homology of Gaussian Random Fields 

In this section we will give the connection between the Euler integral of a func- 
tion and its persistent homology. This will allow us to interpret our computation 
of the expected Euler integral for Gaussian random fields as a computation on 
the expected value of a quantitative measure of a Gaussian random field's persis- 
tent homology. We will start off by giving a brief sketch of persistent homology. 
For more details and further references see [6, 7, 8]. 

6.1. Persistent Homology and the Euler Characteristic 

Persistent homology is a way of tracking how the homology of a sequence of 
spaces changes. For simplicity, in what follows all homology will be with rational 
coefficients. Given a filtration of spaces X = {X u } u such that X s C X t if 
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Fig 1. The persistent homology of a function f : R — > K. In this case the sublevel sets are 
intervals, and therefore the only nonzero homology is Ho, representing connected components. 
The figure to the right presents the barcode of the zeroth persistent homology. A new Ho 
element in the persistent homology is born at each local minimum, and a death occurs at each 
local maximum. When two elements merge, the last one to be born is the first one to die. 



s < t, the persistent homology of X, PH*(X), consists of families of homology 
classes that 'persist' through time. Explicitly an element of PHk(X) is a family 
of homology classes a = {at} for t £ [a,b], where at € Hk{X t ) (the fc-th 
homology group of X t ). These elements are related by the fact that the map 
Hk(X s ) — > Hk(X t ), induced by the inclusion X s C X t , maps a s to a t . The birth 
time a of the element a can be thought of as the first time a appears, which is 
defined by the condition that a a is not in the image of Hk(X s ) — > Hk(X a ) for 
all s < a. The death time b of the element a is the moment that a* becomes 
equivalent to something that existed before a. Formally we require that at is 
not in the image of Hk(X s ) — > Hk(X t ) for all s < a and t < b, but a^ is in the 
image of Hk{X s ) — > Hk(Xb) for all s < a. One must put 'tameness' conditions 
on the filtration so that the birth and death times are defined. 

Given a tame function / : X — > K, there is an associated filtration of spaces 
oo, One defines the persistent homology of / to be the persistent 

homology of this filtration, i.e. PH sr {f) — P_ff > „({/ _1 (— oo, u])}). The persistent 
homology of a tame function / : X M can be seen as a generalization of Morse 
theory, for if / is a Morse function then the critical values will correspond to 
birth and death times of elements in the persistent homology PH*(f). 

A graphical way of representing PH*(X) is via barcodes. Given a persistent 
homology element a, it can be represented by a bar starting at its birth time and 
ending at its death time. To form a barcode, first choose a basis for PH*(X), 
a collection of persistent homology elements such that, for every time t, those 
classes that are alive at time t form a basis for H^(X t ). The barcode for PH*(X) 
will then be the collection of bars for the chosen basis. 

It turns out that the Euler integral of a tame function / is strongly related 
to the persistent homology of /. In light of Proposition 2.4, this is not sur- 
prising, since the Euler integral is a measure of how the Euler characteristic of 
oo, u] changes, while the persistent homology tracks how the homology of 
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oo, u] changes. To make the relationship precise we need to introduce the 
following natural extension of the Euler characteristic to barcodes. 

Definition 6.1. The Euler characteristic of a barcode with a finite number 
of bars and no bars of infinite length is 

bjeB, 

where fx(bj) is the degree ofbj (i.e. the homology degree of the class it represents), 
and £(bj) is the length ofbj. Equivalently, 

x( B *) = / XB(u)du 
Jm 

where Xb{u) — ^2 k (—l) k # {bars of degree k at time u} is the signed sum of the 
number of bars at time u. 

Proposition 6.2. Let f : X — > R be a tame function and let PH*(f, a) be the 
persistent homology of f in the range (— oo,a]. Then 



X (PH*(f, U) = / max x{X)- f\d X ] , 

Jx 

and, in general, 

X (PH*(f,a)) = a X {X)- [ (G a of)\d X l 

Jx 

Proof. Observe that 

/•oo />0 

f\dx\= (x(X)-x(f <u))du~ X{f<u)du 



X JO J-oo 

= /maxXpO - / X(f < U) du. 

J — oo 

Using this equality and \(f ^ u ) — XPH # (f)( u )i by Definition 6.1 we have 
x(PH*(fJ™)) = X(f < u)du = f maxX (X) - [ f\dx\. 

J-oo JX 

As for the second claim, first suppose that a < f max . Then a — (G a o /) max 
and PH*(f, a) = PH lr {G a o f,a), so by applying the first claim to (G a o /) we 
get that 

X (PH*(f,a)) = X (PH*(G a of,a))=a X (X)- f (G a o f)[d X ]. 

Jx 
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If a > / max , then 

X(PF*(/, a)) = [ X (f< u)du = (a - / max )x(X) + f X (f < u)du 

J — oo J — OO 

= a X (X)+ ( f\dx\. 



Jx 

However, a > /max implies that / = G a ° /, so we are done. □ 

6.2. The Expected Euler Characteristic of the Persistent Homology 
of a Gaussian Random Field 

In light of the the connection between the Euler integral of a function and the 
Euler characteristic of the function's persistent homology in place, we will now 
reinterpret our computations about the expected Euler integral of a Gaussian 
random field. This leads to the following result, which as we described in the 
introduction, seems to be the first result giving a precise form for the expected 
value of a quantitative property of the persistent homology of random functions. 

Theorem 6.3. Let f : M —> R fc be a Gaussian random field satisfying the GKF 
conditions, G : R fc — > R continuous and piecewise C 2 , and g = Go/. Then 

E {x(PH* (g, 5max ))} = X (M) (E { 5max } - E {g}) 

d „ 

+ V(2 7 r)- J '/ 2 £ i (M) / Mj{D u )du. 

If f : M —¥ R is a Gaussian random field satisfying the conditions of the GKF, 
then 

E{x(PH*(f, /max))} = E{/ max } X (M) H 1 1 ' - Ul 



V2tt 

Proof. By Proposition 6.2, 

E {x(PH*(g, 3max ))} = E {.g max } X (M) - E U g \d X ] 

Now use the computation of E { J M g \d\~\ } from Theorem 4.1. □ 

One drawback of this result is that it requires knowledge of E{g max }, which 
is usually unavailable. However, for real Gaussian random fields, there is a way 
to circumvent this problem. 

Theorem 6.4. Let f : M — > R be a Gaussian random field satisfying the GKF 
conditions. Then for any a € R, 

E {x(PH*(f, a))} = X (M) (<p(a) + o$(o)) 

d 

+ <p{a) (27r)-- ? '/ 2 £,(M)ff,_ 2 (- a ). 
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Proof. Again, Proposition 6.2 gives that 



E { X (PH. (/, a))} = a X (M) - [ (G a o f) \dx\ 



JM 

and Lemma 5.3 computes J M (G a /) [dx~\ ■ 



□ 



7. An Application 

An interesting application of the Euler integral is suggested in [4]. Suppose 
that an unknown number of targets are located in a space X, and each target 
a is represented by its support U a C X. Suppose also that the space X is 
covered with sensors, reporting only the number of targets each sensor sees (i.e. 
no identification). Let h : X — ¥ Z be the sensor field, i.e. 



The following theorem states how to combine the readings from all the sensors 
and get the exact number of targets. 

Theorem 7.1 ([4]). If all the target supports U a satisfy x(U a ) — w for some 
w 7^ 0, then 



Note that we do not need to assume anything about the targets other than 
they all have the same Euler characteristic. For example, we need not assume 
that they are all convex or even have the same number of connected components. 
On the other hand, the theorem assumes an ideal sensor field, in the sense that 
the entire (generally continuous) space X is covered with extremely accurate 
sensors (the range of each sensor is a single point in X). In [3] more realizable 
models using the lower/upper Euler integralare discussed. 

Using the results from Section 4 we can extend the setup above to the case 
where the readings from the sensors are contaminated by a Gaussian (or Gaus- 
sian related) noise fix). We will use the following proposition. 

Proposition 7.2. let h, f : X — > K be tame functions and suppose that h(X) 
is discrete, then 



Proof. Let h{x) = Y)j—i li^-Aiix), where the Ai are disjoint. Then by the addi- 
tivity of the Euler characteristic we have that 



hix) = # {targets activating the sensor located at x} . 






(7.1) 
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Next, 

/ (h + f)\d X ]= I (a i + f)\d X ]=a i x(A i )+ [ f\d X ], 

JA Z JAi J Ai 

where the last equality follows from Proposition 2.4, since every critical value is 
changed by aj. Applying this to (7.1) completes the proof. □ 

Returning to the target enumeration problem, we have a deterministic signal 
x = j x h\d\\, observed via a noisy measurement Y = j x (h+ f)\d\\. By the 
above proposition we have that 

Y = [ (h + f)[d X ] = [ h\dx\+ I f\d X ]=x + N, 
Jx Jx Jx 

so we have the classical parameter estimation with additive noise model. If 
f{x) is a Gaussian or Gaussian related random field satisfying the conditions in 
Theorem 4.2, then we can use the estimator x = Y—E {N}. Further investigating 
the properties of the Euler integral might lead to useful estimation techniques 
for this model. 
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